Parking in the city: an example of limited resource sharing 
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Abstract 



During the attempt to park a car in the city the drivers have to share limited resources (the available roadside). We 
' show that this fact leads to a predictable distribution of the distances between the cars that depends on the length of the 

street segment used for the collective parking. We demonstrate in addition that the individual parking maneuver is guided 
by generic psychophysical perceptual correlates. Both predictions are compared with the actual parking data collected in 
^ ' the city of Hradec Ki'alove (Czech Republic). 

C/3 ' 

o 1 Introduction 

• I— I ' 

J>-»' Everyone knows that to park a car in the city center is problematic. The amount of the available places is limited and 
, has to be shared between an increasing number of cars. There are several attempts to tackle this challenge by introducing 
parking charges, building underground garages, implementing parking zones, advertising for public transport and many 
others. The problem seems however to persist. Deeper understanding of the related processes is therefore of a common 
interest. 

There have been several mathematical attempts to tackle the parallel parking process. The classical way to do so is 
the "random car parking model" introduced by Renyi (jj - see ||2l, El for review. In this model the cars park on randomly 
, chosen places and once parked the cars do not leave the street. All cars are usually assumed to be of the same length 
T— I ' and the process continues as long as all available free places ale smaller then Iq. The model leads to predictions that can 
, be easily verified. First of all it gives a relation between the mean bumper - to - bumper distance D and the car length: 
D 0.337 ^0- Further: the probabiHty density p{D) of the car distances D behaves like 111,1151, IS, Q p{D) « -ln{D). 
00 , This means that small distances between cars are prefeiTed. 

' To test this results real parking data were collected recently in the center of London [Sl. The average distance between 
J> the parked cars was 152 cm which fits nicely with the relation D ^ 0.337 /q for /q = 450 cm. The predicted probability 
density was however incompatible with the observed facts. The model leads to p{D) oo for D — > whereby the data 
from London display p{D) as D 0. The same behavior has been found also in other cities [|9l- 
■ This results show that the parking process is not so simple as assumed by the "random car parking model". First of 
all: although the total available parking place remains unchanged the cars are reshuffled many times during the day since 
some parked cars leave and new cars park on the vacant places. Moreover the parking maneuver is not trivial. It is not just 
a simple positioning of the car to the parking lot. 

Our aim here is to give and alternative approach to the parking process. We will show that it can be understood as a 
statistical partition of the limited parking space between the competing persons trying to park. The partition is described 
by the Dirichlet distribution with a parameter g. We will also show that this parameter is in fact fixed by the capability of 
the driver to exploit small distances during the parking maneuver 

Let us focus on the spacing distribution (bumper to bumper distances) between cars parked parallel to the curb. We 
will assume that the street segment used for parking starts and ends with some clear and nondisplaceable part unsuitable 
for parking. It can be a driveway or turning to a side street. Otherwise the parking segment is free of any kind of parking 
obstructions. We will assume that it has a length L. Moreover there are not marked parking lots or park meters inside it. 
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So the drivers are free to park the car anywhere in the segment provided they find an empty space to do it. We suppose 
also that all cars have the same length ^o- 

Many cars are cruising for parking in this part of the city. So there are not free parking lots and a car can park only when 
another parked car leaves. To simplify the further formulation of the problem and to avoid troubling with the boundary 
effects we assume that the street segment under consideration form a circle. The car spacing distribution is obtained as a 
steady solution of the repeated car parking and car leaving process. 

Due to the parking maneuver one needs a lot of a length w l.3lo to park. Hence in a segment of length L the 
number of the parked cars equals to iV « [L/(1.3Zo)]- Denoting by Dk the spacing between the car k and fc + 1 we get 
J2k=i = L — NIq and after a simple rescaling finally 

TV 

= (1) 

k=l 

Since all parking lots are occupied the number of parked cars is supposed to be fixed. The repeated car parking and 
car leaving reshuffles however the distances Dk- We will treat them as independent random variables constrained by the 
simplex ([T]). The distance reshuffling goes as follows: In the first step one randomly chosen car leaves the street and the 
two adjoining lots merge into a single one. In the second step a new car parks into this empty space and splits it again into 
two smaller lots. Such fragmentation and coagulation processes were discussed intensively since they apply for instance 
to the computer memory allocation - see [TOl for review. The related equations are simple. If a car leaves the street and 
the neighboring spacings - say the spacings Z3„, Dn+i - merge into a single lot D we get 

D = Dn + Dn+i+lo- (2) 

When a new car parks to D it splits it into Dn,Dn+i- 

bn = aiD-lo) 
bn+i - {l-a){D-la). (3) 

where a G (0, 1) is a random variable with a probability density q{a). The distribution q{a) describes the parking prefer- 
ence of the driver. We assume that all drivers have identical preferences, i.e. identical q{a). (The meaning of the variable a 
is straightforward. For a = the car parks immediately in front of the car delimiting the parking lot from the left without 
leaving any empty space. For a = 1/2 it parks exactly to the center of the lot D and for a = 1 it stops exactly behind the 
car on the right.) Combining (|2]i and (O gives the distance reshuffling 

£>„ = a{Dn+Dn+l) 

bn+1 = (l-a)(A. + (4) 

(The car length drops out.) The simplex ([T]) is of course invariant under this transformation. 

The mappings (HI are regarded as statistically independent for various choices of n. Moreover all cars are equal. So in 
the steady situation the joint distance probability density P{Di, Dm) has to be exchangeable ( i.e. invariant under the 
permutation of variables) and invariant with respect to (|4|i. Its marginals p{D) (the probability densities of the particular 
spacings) are identical: 

Pk{Dk)^p{Dk)= [ P{Du-,DN)dDi..dDk-idDk+i..dDN. (5) 

JDi + ..+Dn = 1 

A standard approach to deal with the simplex (HJ is to regards Dk as independent random variables normalized by a 
sum: 

dk 

Dk = ^rr- (6) 

Here dk are statistically independent and identically distributed and it is preferable to work with them. Moreover: it 
is obvious that the distribution of {Di, ..,Dn} is invariant under the transform (HJi merely when the distribution of 
{di, .., d^} is invariant. So let us apply the relation (HI on the variables c?„: 

dn = a{dn + dn+i) 
dn+1 = {1 ~ a){dn + d„+i) (7) 
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where a, o?„, d^+i are independent and dn, dn+i identically distributed . Giving the distribution q{a) of a we look for 
distributions of dn such that the transformed variables dn, dn+i preserve the distribution of dn - The effort is to solve the 
equation 

d = a{d + d') (8) 

where d' is an independent copy of the variable d and the symbol = means that the left and right sides of (O have identical 
statistical properties. 

Distributional equations of this type are mathematically well studied - see for instance lfT2ll - although not much is 
known about their exact solutions. In particular it is known that for a given distribution q{a) (describing the parking habit) 
the equation ([8]) has an unique solution which can be obtained numerically. Since we are interested in explicit results we 
choose q{a) from a two parametric class of (3 distributions. Then the solution of ([8]) results from the following statement 

ini: 

Statement: Let di, c?2 and a be independent random variables with distributions: di ^ r(ai, 1), ^2 ^ r(a2, 1) and 
a ^ /3(ai, 02). Then a{di + ^2) ^ r(ai, 1). 

(The symbol ~ means that the related random variable has the specified probability density. T{g, 1), (72) denote 
the standard gamma and beta distributions respectively.) 

Since in our case the variables di , ^2 are equally distributed we have gi = 92 = 9 and a ~ (3{g, g). So the probability 
density of a is symmetric in this case, i.e. the variables a and 1 — a have the same distribution. In other words the drivers 
are not biased to park more closely to a car adjacent from the behind or from the front. 

The solution of (O is in this case equal to d ^ r((7, 1). The relation ^ returns the spacings Dk and we find that the 
joint probability density P{Di, Dn) is nothing but a one parameter family of the multivariate Dirichlet distributions 
on the simplex ([T) ifTSll : 

PiD,,...,DN) = ^^Dr'Dr\..D%-' (9) 
Its marginal Q is simply D ~ P{g, {N — 1)5). Normalizing the mean of _D to 1 we are finally left with 

p{D) = ^p(g,{N-l)g,^) . (10) 



Despite of the symmetrical parking maneuver this distribution is asymmetric. This is a consequence of the persistent 
parked car exchange and can be regarded as a collective phenomenon. 

The above considerations leave the parameter g free. But we show that there is in fact a natural choice of g leading to 
g = 3. The point is that the behavior of q{a) for small a reflects the capability of the driver to estimate small distances. 
The collision avoidance during the parking maneuver is guided visually and this ability is shared equally by all drivers. If 
it applies the behavior of q{a) for small a has to be generic, i.e. independent on the particular city or parking situation. It 
is just fixed by the human perception of distance. 

Distance perception is a complex task and there are several cues to do this. Some of them are monocular (linear 
perspective, monocular movement parallax etc.), others oculomotor (accommodation convergence) and finally binocular 
(i.e. based on stereopsis). All of them work simultaneously and are reliable under different conditions - see lfT4ll for more 
details. For the parking maneuver however the crucial information is not the distance itself but the estimated time-to- 
collision between the bumper of the parking car and its neighbors. This time has to be evaluated using the knowledge 
of the distance and velocity. It has been argued in a seminal paper by Lee IfTSl that the estimated time to collision is 
psychophysical^ evaluated using a quantity named r . It is defined as the inverse of the relative rate of expansion of 
the retinal image of the moving object. Behavioral experiments have indicated that r is indeed controlling actions like 
contacting surfaces by flies, birds and mammals (including humans): see 1 16 Hf7l . l|T8l . 

When the observer moves forward in the environment, the image on the retina expands. The rate of the expansion r 
conveys information about the observer's speed and the time to collision. Psychophysical and physiological studies have 
provided abundant evidence that r is processed by specialized neural mechanisms in the brain [19]. We take r to be the 
informative variable for the final braking - see [20 1 for review. 

Let 9 be the instantaneous angular size of the observed object (for instance the front of the car we are backing to 
during the parking maneuver). Then the estimated time to contact is given by 
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Since 9{t) = 2 a.rctaii{LQ /2D{t)) with Lq being the width of the approached object and D{t) its instantaneous 
distance, we get 

'^(^) = -TTT^TTTTTTTTT arctan — ^ . (12) 



2LoidD{t)/dt) \'2D{t)^ 

For D >> Lq and a constant approach speed v — —dD/dt the quantity r simply equals to the physical arrival time: 
r = D/v. For small distances (parking maneuver), however, r « D^/ (wLq) and the estimated time to contact decreases 
quadratically with the distance. 

Let us return to the equation dUi. For a fixed parking lot the final stopping distance is proportional to a. Assuming that 
the courage to exploit small distances is proportional to the estimated time to contact we finally get for the probability 
density q{a): q{a) « for small a. Since q{a) = f3{g, g, a) the behavior q{a) « fix the parameter g to g = 3 and the 
normalized clearance distribution ( fTOl i reads 

piD) = 1/? fa, 3(iV '^"""^ -I^^D^^N - Dr-^ (13) 



iV"^ V \N J 2r(3(7V-l)) 

so the parameter g is fixed to 3. But the distribution still depend on the number of cars N in the parking segment. 
It is a consequence of the constrain ([T]i. For large number of cars, >> 1, the constrain ([T]i does not play a substantial 
role and p{D) equals to r(3, 1, D) (this is true in the limit N — > oo). In the other extreme case with = 1 (the parking 
segment is so short that it allows the parking of a single vehicle) the distribution p{D) just reflects the parking maneuver 
and is equal to q{D). For parking segments of intermediate size the theory predicts a dependence of the results on the 
segment length. 

To verify the predictions of the model we measured the bumper to bumper distances between cars parked on two 
different streets in the center of Hradec Kralove (Czech Republic). Both streets were located in a place with large parking 
demand and usually without any free parking lots. In addition one of these street (street 1) contained driveways to court- 
yards. This means that the actually available fixed parking segments were much shorter on this street. In the mean 3-4 cars 
were able to park among two subsequent driveways and we collected 773 car spacings under this conditions. The second 
street (street 2) was free of any dividing elements. Here we measured altogether 699 spacings. 

The probability distributions resulting from these data seems to be fairly compatible with the prediction of the model. 
First of all: the perceptual mechanism based on the estimated time to contact seems to be verified. We demonstrated 
that when the estimated time to contact is decisive for the final car stopping then the parameter g in ( fTOb fix to g ~ 3. 
And exactly this value fits with the measured data. Moreover: the finite length of the street (the simplex ([T])) leads to a 
dependence of the spacing distribution (T3[ on A^ . So the result obtained for short and long parking segments should be 
different. And this is indeed observed when the data from the street 1 and 2 are compared. We plot the results on the figure 

m 

The difference between the results is not large but it is nevertheless clearly visible (compare the full and dashed lines). 

To summarize we have shown that the clearance distribution for the cars parked in paralel can be described as a 
marginal distribution of the multivariate Dirichlet distribution with a parameter g. The parameter is fixed to g = 3 by the 
psychophysical^ estimated time to collision during the parking maneuver. The measured data support this hypothesis. 
The theory leads further to a prediction that the clearance distribution depends on the length of the used parking segment. 
Also this fact is verified by the collected data. 
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Figure 1 ; The measured distance distributions for cars parking on the the street 1 and 2 in the city of Hradec Kralove are 
compared with the prediction of the formula il3[ . The resuhs for the street 1 and 2 are plotted on the panels a) and b) 
respectively. Bars mark the probability density extracted from the collected data. The full line stays for the formula ( fTsT i 
with iV = 3 (panel a)) and N=20 (panel b)). To guide the eyes the result of (fTsT i for = 20 and = 3 are plotted on the 
panels a) and b) as a dashed line. 
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